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We study the nature of the ground state of the quantum 
dimer model proposed by Rokhsar and Kivelson by diago- 
naUzing the Hamiltonian of the model on square lattices of 
size L X L, where L < 8, with periodic boundary conditions. 
Finite-size scaling studies of the columnar order parameter 
and the low lying excitation spectrum show no evidence of a 
dimer liquid state in any finite region of the zero temperature 
phase diagram. In addition, we find evidence of a transition 
from the columnar dimer state to an intermediate state at a 
negative value of V/ J. This state is identified to be the pla- 
quette resonating valence bond (RVB) state. The energy gap 
of this state vanishes as a power law of L. It exhibits colum- 
nar dimer order, but has disorder within the dimer columns. 
This state persists up to V/ J < 1, and the system changes to 
a dimer liquid state only ai V/ J = 1. 

PACS: 75.10Jm, 75.40Mg 



I. INTRODUCTION 

The quantum dimer model (QDM) was first introduced 
by Rokhsar and Kivelson (RK)EJ. It was proposed as an 
alternative description of the non-Neel state of the spin-i 
Heisenberg antiferromagnet on a square lattice. In a state 
with exponentially decaying spin-spin correlation and a 
large energy gap separating the spin excitations from the 
singlet ground state, RK argued that the low energy 
physics is contained in the "short-range resonating va- 
lence bond (RVB)" states, spanned by the set of nearest- 
neighbor valence bond states. Nearest-neighbor valence 
bonds are called dimers. The basis set thus consists of 
all possible dimer configurations at closest-packing (all 
sites form exactly one dimer with one of their nearest 
neighbors). RK also argued that the non-orthogonality 
of these basis states can be absorbed by defining a phe- 
nomenological Hamiltonian with short range dimer inter- 
actions. The relevance of this quantum dimer model to 
the frustrated quantum antiferromagnet was supported 
by other injdependent studies. Large- expajision,p series 
expansion,0 and numerical diagonalizationP have shown 
that a spin-Peierls state exists in some region of the phase 
diagram where the Neel state is unstable. This spin- 
Peierls state corresponds to the state in the QDM where 



the dimers are frozen into a columnar pattern (Fig. 0(a)). 
Recently the QDM has been generalized to the kagome 
lattice by Zeng and Elser.B 

The Hamiltonian of the QDM proposed by RK is 
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II and — represent parallel dimers on the opposite sides 
of a square plaquette. The first term of Ti. is the dimer 
kinetic energy operator while the second term is the po- 
tential energy. The space spanned by all close-packed 
dimer configurations on a square lattice can be divided 
into distinct topological sections characterized by a pair 
of conserved winding numbers (QxT^y) (refer to Ap- 
pendix^. On an Lx L square lattice, the allowed wind- 
ing numbers are —L/2 < i^x,^y < L/2. At V/J > 1, 
the exact ground state of the QDM is the staggered 
dimer solid shown in Fig. 0(b). It is four- fold degenerate 
and has winding numbers (±L/2,0) and (0, ±L/2). At 
V/J ^ —1, the ground state develops columnar order as 
shown in Fig. |l|(a). At V/ J = 1, the lowest energy states 
in all topological sectors are zero-energy eigenstates of 
Ti.. They are equal superpositions of all the dimer config- 
urations in their sectors. Therefore, the ground state 
ai V/J = 1 is a dimer liquid, and any ground state 
dimer correlation function can be calculated exactly by 
the method of Fisher and StephensonO for the classical 
dimer problem. Consequently this state is called the FS 
state. The ground state properties of the QDM in the 
range < V/ J < 1 are not so clear. There are three 
main questions: (a) is the ground state in this range or- 
dered (dimer solid) or disordered (dimer liquid)? (b) If 
it is ordered, what is the nature of this order? And (c) if 
it is disordered, is there a gap in the low lying excitation 
spectrum? RK suggested that a dimer liquid state can 
exist in a finite region ofV/ J < 1. By numerically diagctj 
nalizing the QDM on square lattices up to 6 x 6, SachdevCl 
found evidence for the columnar state, but no evidence 
for a dimer liquid state at V/J ^ 1. However, due to 
finite-size effects, he did not rule out the possibility that 
dimer liquid exists over a finite range k < V/ J < 1, and 
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estimated a lower bound 0.5 < k. On the other hand, 
some analytic studies have suggested the existence of a 
dimer liquid state at V/ J ^ I. By mapping the QDM 
to a roughening problem, LevitovB showed that|-a dimer 
liquid with gapless excitations may exist. OrlandQ formu- 
lated the QDM as a system of non-interacting fermionic 
strings, and showed that at V/J — 0, the ground state 
is a dimer liquid. From this he inferred that the ground 
state is a dimer liquid for < 1^/ J < 1, with a phase tran- 
sition to a dimer sohd at some V/J < 0. We also note 
that another kind of dimer order has been proposed,t2l 
although in a different context. This is a plaquettc RVB 
state as shown in Fig. |l](c) , where every other plaquettc 
is in the | || ) or | =) state with equal probability and in- 
dependent of other plaquettes. In this paper, we extend 
Sachdev's calculations to an 8 x 8 square lattice. We aim 
to identify the different phases of the QDM in the range 
V/J < 1. 

II. NUMERICAL CALCULATIONS 

Using the Lanczos algorithm, we diagonalize the 
Hamiltonian H in the range — l<T^/J<lonanLxL 
square lattice with periodic boundary conditions, where 
L = 4, 6, and 8. The basis states are all possible dimer 
configurations at closest packing. TJae number of basis 
states for a given L can be evaluatedQ analytically. Enu- 
merating all the dimer configurations is much more diffi- 
cult. Nevertheless, we have devised an efficient algorithm 
which enumerates all the 300 million dimer configurations 
on the L = 8 lattice within four hours using an HP 735 
workstation (see Appendix ^) . Using translational sym- 
metry, we can reduce the number of basis states on an 
L X L lattice by a factor of about L^. For V/J < 1, 
the ground state has momentum k = (0,0), and wind- 
ing numbers (0,0). By restricting the basis states to the 
topological sector with winding numbers (0, 0), the num- 
ber of basis states can be further reduced by a factor 
of about 2. To calculate the ground state eigenvector, 
we restrict ourselves to the subspace with momentum 
k = (0,0) and winding numbers (0,0). The number of 
basis states in this subspace is about 2.4 million. This 
makes the numerical diagonalization possible on a work- 
station. 



III. SYMMETRIES OF LOW LYING EXCITED 
STATES 

To study the different phases in the range — 1 < V/ J < 
1, we first study the symmetries of the low lying excited 
states. They are important in finite-size studies. If a sys- 
tem possesses a broken symmetry in the thermodynamic 
limit, the ground state of the finite system will still be 
totally symmetric. In this case the ground state expec- 
tation of the appropriate order parameter will have long- 



range correlations, and there will exist low lying excited 
states with the appropriate symmetries. The columnar 
dimer state shown in Fig. |l|(a) has winding numbers (0, 0) 
and is four-fold degenerate. These degenerate states can 
be combined to form four states, two with momentum 
(0,0), and two with momenta (tt, 0) and (0,7r). Conse- 
quently in a finite system which possesses columnar order 
in the thermodynamic limit, the state with momentum 
(0, 0) will be the ground state and the others will appear 
as low lying excited states, which are degenerate with 
the ground state in the thermodynamic limit. Note that 
the same is true for the plaquettc RVB state shown in 
Fig. |I|(c). Therefore, low lying excited states with mo- 
menta (tt, 0) and (0, tt) may indicate that columnar or 
plaquettc RVB order exists in the ground state. 

We calculate the energies of the lowest few eigen- 
states of the QDM on the L — 4, 6 and 8 lattices for 
— 1 < V/J < 1. All three lattices show qualitatively 
the same picture. The ground state Eq always has zero 
momentum and winding numbers, as mentioned in sec- 
tion ||. Except when V/J is close to 1, the first ex- 
cited state is degenerate with momenta (7r,0) and (0,7r), 
and has zero winding numbers. We call this state q)- 
In the same region, the next excited state is degenerate 
with momenta (0, 0) and (tt, tt), and has winding numbers 
(±1,0) and (0, ±1). We call this state E(^^^^y In Fig. | 
we plot the energies of the states Eq, -E(^,o) ^^id -E'(Tr,7r) 
in a small range of V/J close to 1 for the L = 8 system. 
We can see that there exist a. Sl (which depends on the 
size L) close to but less than 1, such that at V/J = 6l, 
tt) Bjnd £'(7r,o) cross each other and £'(7r,7r) becomes 
the first excited state ioi 6l < V/ J < 1. This crossing 
of energy levels may imply the existence of a new ground 
state in the region 6l < V/ J < 1. However, finite-size 
extrapolation analysis on Sl shows that this is not the 
case._As shown in Fig. 0, Sl approaches 1 as a power law 

ofi,Er 

(1-(5l)(xL-213. (2) 

In other words, such a spurious energy level crossing is 
only a finite-size effect. It does not occur in the thermo- 
dynamic limit and no transition is likely to occur at V/J 
close to 1. Therefore, we expect that the -E(7r,o) state is al- 
ways the lowest excited state in the range —1 < V / J < 1. 
This is consistent with, but does not imply, the existence 
of long-range order (LRO) in the QDM at -1 < V/J < 1. 

Next we study the energy gap AE^, which is the en- 
ergy difference between -E(^,o) ^nd Eq. Fig. ^ shows dif- 
ferent plots of AEl vs L at different V/J. When V/J 
is close to 1, complication arises because of the spurious 
crossing over of the states £'(,r,o) and -E(^.Tr) a.t Sl- But 
we find that AEl is not very sensitive to V/ J in the re- 
gion close to 6l (see Fig. ||). Consequently we plot AEl 
at V/J = 6l (instead of at the same V/J for different 
L) when V/J is close to 1. From Fig. Q we can clearly 
distinguish two different behaviors. At V/J = —1, AEl 
decays exponentially with L. This is a clear signal for 
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the existence of LRO. When < T^/J < 1, we find 
that Ai?i vanishes with a power law in L. Best fit to 
the data shows that the power is approximately 2 in the 
whole range. For —1 < V/J < 0, we are not able to 
determine whether AE^ has exponential or power law 
dependence on L. This is probably because of the more 
serious finite-size effect when V/J is close to a transition 
point. Hence we infer that there are two phases in the 
range —1<V/J<1. At some negative V/J, a tran- 
sition occurs and the dependence of AEl on L changes 
from exponential to power law. 



G(r) is a nonzero constant at large r. When V/J is close 
to 1, it is not possible to tell from Fig. ^ whether G{r) is 
small or zero at large r. Nevertheless, it tells us that in 
the range ~1 < V/ J < 1 the long-range order, if exists, 
should be columnar. 

As suggested in Ref. a further probe of the columnar 
dimer order is provided by the cumulant of the columnar 
order parameter, 
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IV. COLUMNAR ORDER PARAMETER 



The columnar dimer order parameter is defined asE 
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where x and y are unit vectors. The dimer number oper- 
ator n(r-|-|) is 1 if the site at r and its nearest neighbor at 
r -|- e form a dimer, and zero otherwise|^ In the finite-size 
study, one defines the order parameteifl 
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For a dimer liquid which has no long-range dimer order, 
XL is zero in the large L limit. If long-range columnar 
order exists, xl remains finite at large L. In Fig. ^ we 
plot XL versus V/J. When V/J is close to —1, xl is ap- 
proximately linear in 1/ L^. But when V/J is increased, a 
linear relation in 1/L fits the data better. Unfortunately 
we are not able to determine the point where these two 
behaviors change from one to another. Using these linear 
relations, we extrapolate xl to obtain Xoo, which is also 
plotted in Fig. ||. It shows that Xoo is finite at V/J < 1, 
which signifies the existence of columnar order. Xoo de- 
creases as V/ J approaches 1 , showing that the columnar 
order is weakening. But there is no evidence that Xoo is 
zero at any V/J < 1. Xoo clearly shows that the columnar 
order persists for V^/J < 1. 

Different dependence of xl on L (oc 1/L^ and 1/L) 
may imply different orders in the respective regions. One 
defines the correlation function of 5'co/|3 
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where r, ri are in the A sublattice. In the two regions 
where xl has finite-size corrections and 1/L, G{r) 

probably decays exponentially and as a power law in r, 
respectively, to a non-zero constant value. Fig. ^ shows 
G{r) in the L = 8 system. Unfortunately, direct obser- 
vation of how G(r) decays with r is not possible in this 
system size. Fig. ^ shows that when V/J is away from 1, 



If long-range columnar dimer order exists, — *■ 1 as 
L ^ 00. We plot gL vsfsus V/J in Fig. 0(a). Standard 
finite-size scaling theorylla shows that at phase transition, 
(7i at different L cross at a unique point. Since gL do 
not cross except when V/J is close to 1, Fig. ^a) tells 
us that the order parameter '^coi detects no transition 
except when V/J is close to 1. This is the transition to 
the dimer liquid (FS) state. Since our system size L is 
not large enough to make cross at one point, it is not 
possible to determine the transition point from Fig. ^a) . 
Nevertheless, our previous finite-size analysis oi AEl and 
Xl strongly suggest that some kind of columnar LRO 
persists all the way up to T^/ J = 1. We also note that g^ 
detects no transition at any —1 < V/J < 0, in contrary 
to the results of AEl and xl- 



V. DISORDER WITHIN THE COLUMNS 

It is easy to understand why gj^ is not able to detect 
the transition at V/J < 0. Note that in the plaquette 
RVB state, Xi is also non-zero in the large L limit. It 
can be considered as a state which possesses long-range 
"columnar order" (as measured by xl), but has disorder 
within the established columns. To detect this disorder 
we introduce an order parameter Mm ^, 
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where n{\ \\ )) and n{\ —)) are number operators of ver- 
tical (II) and horizontal {—) dimer pairs respectively. In 
the perfect columnar state, |Af||^| = 1/2, while in the 
plaquette RVB, dimer liquid and the staggered state, 
-^^11.= = 0. Analogous to we define the cumulant 
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Again, finite-size scaling theory shows that at the phase 
transition, g\\^^ at different L should intersect at a unique 
point. In Fig. ^(b) we plot g\\- vs V/J for different L. 
Although the three curves do not intersect at a single 
point, the crossing of the curves at L = 6 and 8 indicates 
that a transition is possible at some negative value of 
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VjJ. This is consistent with the results of A£^i and 
XL- We are not able to obtain a good estimation for the 
transition point Vd J ■ But very roughly, Vd J ^ —0.2. 
We also calculate the moments of M\ _ , 
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i.e. the difference between the number of vertical and 
horizontal dimers. The corresponding cuniulant g\^_ is 
shown in Fig. 1(c). We find a similar crossing of g\- 

for i = 6,8 also near VjJ 0.2. Both M\_ and M\\- 

should be probes of the same ordered columnar to plaque- 
tte RVB transition. We note that for V > Vc both 
and 51 are close to 3. This value is the expected gaus- 
sian moment ratio for disordered states. This strongly 
suggests there is no or M| LRO for the whole 

range oiV>Vc. For V/ J < —1 the cumulants approach 
1 for the ordered columnar state. In the L —>■ oo limit 
the cumulants should become 1 for all V <Vc- Fig. ^(b) 
and (c) are consistent with this scenario. 



VI. DIMER CORRELATIONS 

311^= and .g| _ show that the possible transition at Vc/ J 
is between two states with long-range columnar order, 
but one has disorder within the columns. Next we are 
going to study the dimer correlations in these two states. 
The spatial correlation of the dimers can be displaj^ed by 
evaluating the dimer-dimcr correlation function,E3 
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where = n ( ^ (r^ +r ^ ) ) . If the dimers do not have LRO , 
C(ij)(ki) should show some short-range correlations, and 
fall off rapidly as the separation between the bonds {ij) 
and [kl) increases. On the other hand, if the dimers have 
LRO, C(ij)(fe;) should reflect the pattern of the long-range 
correlations. Fig. ^ shows the dimer-dimer correlation in 
the L = 8 system at V/ J = —1 and 0. The reference 
bond {ij) is represented by a double line. Ci^ij-^(^i^i^ is pro- 
portional to the thickness of the line joining the pair of 
sites k and I. Solid line means C(ij)(M) > (correlation), 
and broken line means C(^ij)(^ki) < (anti-correlation). It 
is obvious that in both cases the overall order is colum- 
nar — those bonds with positive ^(^^(fc;) are arranged in 
well-defined columns. This is consistent with the study 
using 5'co/(r) as the order parameter (Sec. Except 
for the trivial short-range correlations, C(^ij)(^j^ij does not 
fall off significantly with distance. This is an indication 
that the correlation is long-range. A major difference 
between Fig. ||(a) and (b) is that in the latter, C(ij)(fc;) 
for vertical (kl) are very close to zero, except for those 
in contact with (ij). In the perfect columnar state as 
shown in Fig. [^(a), C(^ij)(ki) = 1 or —1/3 when the bond 



{kl) forms and does not form dimer respectively. (Here 
the reference bond (ij) has been chosen to be one of the 
horizontal bonds that form dimer in Fig. |l|(a).) This is 
to be compared to the result dX V/J = —1 (Fig. p|(a)). 
In the perfect plaquette RVB state as shown in Fig.ll^(c), 
except for those {kl) which are in the same plaquette as 



for vertical {kl), and C| 



(y)(fei) 



±1/3 



for horizontal {kl), depending on whether {kl) belongs 
to any one of the plaquettes shown in Fig. |^(c). In the 
same plaquette as {ij), horizontal and vertical {kl) have 
(^{ij)(ki) =2/3 and —1/3 respectively. This is to be com- 
pared to the result at V/J = (Fig. H(b)). Hence the 
dimer correlations strongly suggest that at V/J = —1, 
the ground state is ordered columns, whereas a,tV/J = 0, 
it is the plaquette RVB state. 



VII. PLAQUETTE RVB STATE AND THE 
TRANSITION AT K/J 

From the above results and discussions, there is evi- 
dence that at Vc/ J, the ground state of the QDM changes 
from the ordered columnar to the plaquette RVB state. 
Since this transition does not appear to have been dis- 
cussed in the literature, we will discuss it in a little more 
detail. One can motivate the existence of this transition 
by a simple argument where the QDM is approximately 
mapped onto a two-dimensional Ising model in a trans- 
verse field. The T = transition from ferromagnet to 
paramagnet in that model corresponds to the Vc/ J < 
transition of the QDM discussed in the previous sections. 
We now outline the steps and implications of this approx- 
imate treatment. 

First consider the case V/J <^ —1 where the system 
tends to the ordered columnar state in Fig. |l](a) . The ba- 
sic quantum fluctuations in the ground state are induced 
by the kinetic energy term in (|l|), namely by flipping I =) 
to I II ). For each dimer populated column in Fig. |l|(a) 
consider every other plaquette (as in Fig. ^c)) to have 
an index "i" and state variable af — 1 ii it is in the state 
=} and erf = -1 ff it is | || }. Fig. 0(a) has all af = 1, 
i.e. a ferromagnetic arrangement, whereas Fig. |l|(c) has 
all erf = ±1 with probability 1/2, i.e. a paramagnetic ar- 
rangement. Note also that the other ferromagnetic Ising 
configuration, with all erf = —1, gives another perfectly 
ordered columnar state. The collection of states gener- 
ated by all possible {erf} is, of course, a small subset of 
the set of all close-packed dimer configurations. How- 
ever, they at least describe the plaquette RVB states, 
the ordered columnar states, and one half of the basic 
hipping excitations. Given that our earlier numerical re- 
sults suggest columnar order for all V/J < 1, we do not 
worry too much that the {erf} subset of dimer states are, 
by construction, constrained to have long-range colum- 
nar order, since we are concerned with the disordering of 
other degrees of freedom {e.g. dimer pair arrangements 
and correlations within each established column). 
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Limiting the system to the collection of all possible 
{erf} values, one can show that the Hamiltonian of the 
QDM is, to within an additive constant. 



(11) 



where the sums are over every other plaquette as shown 
in Fig. |l|(c). erf and af are Pauli operators. The kinetic 
energy term with J now induces flipping of the vari- 
ables (like a transverse Ising field hx), and the potential 
energy term V introduces interactions (like an Ising ex- 
change integral iTzj^ between neighboring plaquettes as 
shown in Fig. |l|(c).li3 

For V/ J ^ —1 all the erf are arranged ferromagneti- 
cally, whereas for some K/J < the system undergoes 
a continuous phase transition (in the universality class 
of the classical 3D Ising model) to a state with no LRO 
in the af variables. Indeed, the case F = is trivially 
seen to be perfectly disordered with respect to the erf 
variables, and is precisely the state represented schemat- 
ically in Fig. |l](c). For Vc < V < Q there is some short 
ranged er^ order, but no LRO. Since our (=) (||) pairs of 
dimers are neither allowed to break apart nor to move 
to neighboring plaquettes, this approximation has long- 
range columnar order (as measured by '^coi) preserved 
past the transition point Vc- Our numerical data sug- 
gests this to be the case for the QDM as well: LRO in 
'^coi exists for all V/ J < \ but there is a loss of LRO in 
M\\- at some negative value of V . 

One may speculate that the V/ J <Q transition in the 
QDM is also in the same universality class of the 3D 
Ising model. Unfortunately the system sizes we can deal 
with are too small to determine Vc accurately, let alone 
to estimate the critical exponents. 



VIII. FINITE TEMPERATURE PHASE 
DIAGRAM 

Although our numerical calculations are limited to the 
T = ground state of the QDM, it is interesting to spec- 
ulate on the behavior of the system at finite tempera- 
tures. The simplest, and we feel most likely, scenario 
is displayed schematically in Fig. ^. For any fixed ra- 
tio 1^/ J, as T — » cx) the only interaction that matters 
is the hard-core dimer constraint, and so the FS dimer 
liquid will always be obtained at high enough tempera- 
ture. We expect that for fixed V <Vc as T \s increased 
the ordered columnar state will give way to a plaque- 
tte RVB state,t3 and at a second higher temperature the 
columnar order will be destroyed as the system enters 
the high-temperature dimer liquid phase. Similarly, for 
V > J the staggered phase in Fig. ^ (b) will be destroyed 
at high enough T. The simplest guess is that it goes di- 
rectly to the dimer liquid,E£l although it is not possible to 
rule out other phases. It seems likely that the high tem- 
perature dimer liquid region makes it all the way down 



to the T = axis, but only at the point V = J. It is 
intriguing that the T = oo Fisher-Stephenson (FS) state 
is exactly regained at T = at a single point. 



IX. CONCLUSION 

To conclude, we find no evidence of a dimer liquid state 
in the QDM in the range —1<V/J<1. In this range, 
the QDM exists in two states. At -1 < V/J < Yd J , 
the ground state possesses long-range columnar dimer 
order as suggested in Fig. |l|(a). Its columnar order cor- 
relation function (5'co/(0)5'co/(r)) is likely to decay to a 
constant value exponentially with r, and it has expo- 
nentially vanishing energy gap. This state should exist 
at any VjJ < —1, because negative V favors parallel 
dimers. At Vc/J < V/J < 1, the ground state is the pla- 
quette RVB state as illustrated in Fig.|l](c). Its columnar 
order correlation (5'co/(0)^'co;(r)) is likely to decay to a 
constant value as a power law of r, and the energy gap 
vanishes as a power of L. By approximately mapping 
the QDM to an Ising model with a transverse field, we 
infer that the transition at Vc/J is continuous and pos- 
sibly in the same universality class as the 3D classical 
Ising model. Unfortunately, our system sizes do not al- 
low us to determine the precise value of T4/ J. Our rough 
estimate is Vc/J ~ —0.2. We note that we cannot rule 
out additional transitions in the range Vc/J < V/J < 1, 
we can only say the quantities we have measured do not 
indicate any more. |-, 

Compared to earlier workQ which concluded that a 
dimer liquid state may exist in k <V/ J < 1 and gave an 
upper bound of 0.5 to k, we are able to perform more re- 
liable finite-size scaling analysis with the addition of the 
L = 8 results. Our results push k to very close to 1 which 
suggest that k is in fact 1. Note that although our results 
disagree with Ref. |, our conclusions have something in 
common. Ref. ^ claims that the QDM is a dimer liquid 
a.t V/ J = 0, and exhibits a transition from the columnar 
to the dimer liquid state at some negatipie V/J. We also 
find a transition at some negative j£j However, this 
transition is not from the columnar to dimer liquid state, 
but rather to an intermediate state which we identify to 
be the plaquette RVB state. 
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APPENDIX A: WINDING NUMBERS 

Following Ref. ^ we use the columnar state as the 
reference configuration to define the winding numbers. 
In Fig. |l^(a), the dimers of such a reference configura- 
tion are indicated by dash lines. The arrows along the 
dimers always point from the A sublattice (solid circles) 
to the B sublattice (open circles). In Fig. |l^(b), a dimer 
configuration (indicated by solid lines) is superimposed 
on the reference configuration, but the arrows along the 
dimers (solid lines) point from the B sublattice to the 
A sublattice. This results in directed close paths which 
may wrap around the boundaries (with periodic bound- 
ary conditions). The winding numbers Qx and Qy are 
the net number of loops (clockwise minus counterclock- 
wise) wrapping around the boundaries in the x and y 
directions respectively. Each dimer configuration has a 
unique pair of winding numbers. By repeatedly flipping 
parallel dimer pairs (i.e., applying Ti.) of a dimer configu- 
ration, one can generate all dimer configurations with the 
same winding numbers. But no two configurations with 
different winding numbers are related by such dimer fiip- 
ping operations. Therefore, the Hilbert space spanned 
by all close-packed dimers can be divided into subspaces 
labelled by winding numbers. 

To calculate the winding numbers of a given dimer con- 
figuration, we start from a site and traverse the lattice 
following the arrows as described above until we return 
to the starting site. We mark all the sites along this 
path and count the net number of times it crosses the 
boundaries in the x and y directions. This procedure is 
repeated starting from an unmarked site until all sites are 
marked. It is easy to see that the columnar dimer state in 
Fig. 0(a) has winding numbers (0,0), and the staggered 
state in Fig. [|(b) has winding numbers {zLL/2, 0). 

APPENDIX B: ENUMERATING DIMER 
CONFIGURATIONS 

In this Appendix we describe the way we generate all 
the close-packed dimer configurations on an i x L lattice. 
Again we divide the square lattice into two sublattices A 
and B. At close-pack, each site is one end of exactly one 
dimer. Consequently, a dimer configuration can be rep- 
resented by specifying the direction (±x or ±y) of the 
dimer attached to each site in the A sublattice. There 
are 2 such combinations. But most of them are for- 
bidden because some sites in the B sublattice have zero 
or more than one dimer attached. The most straightfor- 
ward way to enumerate all the dimer configurations is 
to scan through these 2^ combinations to identify the 
allowed configurations. But this takes prohibitively long 
for L = 8. We need a clever way to identify and eliminate 
the forbidden configurations quickly. 

We further divide the A sublattice into two quarter- 
lattices ^0 and Ai. Each site takes up two bits in a 



computer word, with the four binary values 00, 01, 10, 
and 11 representing the direction of the dimer attached to 
it. We divide our task into two steps. First we generate 
allowed configurations for alone, without considering 
Ai. Then we generate the allowed configurations for Ai 
under the constraints of the Aq configurations. The Aq 
and Ai configurations together give all the allowed dimer 
configurations. 

^0 has sites. Each Aq configuration is repre- 

sented in the computer by an integer word with at least 
2 X (LV4) = bits. The 2jth and 2j + 1th bits 

of the integer word represent the direction of the dimer 
attached to the jth site of • In this way each Aq config- 
uration is uniquely associated with an integer value. To 
scan through all the possible configurations, we start 
from an integer value 2^ — 1 whose binary representa- 
tion has Is in the lower L'^/2 bits, and Os otherwise. This 
integer value corresponds to the Aq configuration with all 
dimers pointing in the same direction. Successively de- 
creasing this integer value by 1 is equivalent to rotating 
the dimers in all possible directions starting from the 0th, 
1st, ... sites. In this way all possible configurations 
can be generated. For each integer value (i.e., possible 
Aq configuration), we check for conflict starting from the 
last site. If at some stage we find that the dimers at- 
tached to the zth and jth sites {i > j) of are confiict- 
ing (both connect to the same site in the B sublattice), 
we can skip all the bits lower than the 2jth bit, i.e., in- 
stead of decrementing the integer value by 1, we can go 
to the next possible configuration by changing the 2jth 
and 2j + 1th bits. This substantially decreases the run 
time. For L — 8, Aq has 9 983 558 allowed configurations. 

The set of allowed configurations for Ai is the same as 
that for Aq. Naively we can match these two sets of con- 
figurations and eliminate the forbidden ones. But this 
involves (9 983 558)^ - lO^"* steps for L = 8 and will take 
prohibitively long. We must fully utilize the geometri- 
cal constraints imposed by the configurations when 
enumerating the possible Ai configurations. For a par- 
ticular Aq configuration, we mark the allowed dimer di- 
rections for each site in Ai under the constraints imposed 
by that Aq configuration. The total number of combina- 
tions is substantially smaller than 2^ We then scan 
through all these combinations in the same manner as in 
enumerating the Aq configurations. Repeating the same 
procedure generates all the allowed dimer configurations. 
For L = 4, 6, and 8, the number of allowed dimer con- 
figurations are 272, 90176, and 311853 312 respectively. 
For L = 8, all the allowed dimer configurations can be 
enumerated in about four hours using an HP 735 work- 
station. 



* Electronic address: phleung@usthk.ust.hk 
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^* We note that if the "expansion" is performed about stag- 
gered columns (every other dimer column shifted by one 
plaquette) additional terms proportional to V'^al will 
appear in Ti, but this does not change our conclusions that 
a disordering transition occurs for some V/ J <0. 
Recalling our simple transverse Ising model description, 
this T > transition to plaquette RVB state should be 
in the 2D classical Ising model, not 3D. Even if the trans- 
verse Ising picture breaks down, it should be generally true 
that the quantum phase transition of the QDM at T = 0, 
1/ = Vc is in a 2+1 dimensional classical system universal- 
ity class, while the rest of the T > critical line is in some 
2-dimensional system universality class. 
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Indeed, our analysis cannot rule out 14 = 0. 

FIG. 1. The (a) columnar, (b) staggered, and (c) plaquette 
RVB states. 

FIG. 2. Energy levels _Bo, -E(^,o)i and -B(,r,ir) for the L = 8 
system. All other energy levels are higher than these three 
and are not shown. 

FIG. 3. (1 - 5l) vs L. The straight line is the best fit to 
the data, 1-5l= 3.1544 L'^ ". 



FIG. 4. (a) AEl vs L in semi-logarithmic scale. The 
straight line is the best fit to the data, /^El = 2.742 e"°-^"^. 
(b) Similar plot but in logarithmic scale. 8l is close to 1 (see 
Fig. The straight lines are best fits to the data of the form 
AEl oc1/L^ 



FIG. 5. XL vs V/J for L = 4, 6, 8. A and v are Xoo 
obtained by extrapolating xl linearly in and 1/L re- 

spectively. 



FIG. 6. G(r) at different V/J for L = 8. 

FIG. 7. Various cumulants vs V/J for L = 4, 6, and 8. (a) 
gL, (b) -, and (c) g\_. 

FIG. 8. Dimer-dimer correlation C\ij^ki) in the L — 8 sys- 
tem at (a) V/J = —1 and (b) V/J = 0. Only a quadrant 
of the system containing all the inequivalent dimer pairs is 
shown. The reference bond (ij) is represented by a double 
line. Ci^ij){ki) is proportional to the thickness of the line join- 
ing the pair of sites k and I. Solid line means C{^ij)(ki) > 0, 
and broken line means C(^ij)(M) < 0. Note that the line widths 
are in the same scale in both cases. The values of C{ij)(u) are 
shown next to the bonds, where different fonts are used for 
vertical and horizontal bonds. 



FIG. 9. Speculated phase diagram of the QDM. 

FIG. 10. (a) The reference columnar dimer configuration, 
with arrows pointing along the dimers from the A sublattice 
(•) to the B sublattice (o). (b) A dimer configuration (in- 
dicated by solid lines) superimposed on (a), but with arrows 
pointing in the opposite direction, i.e., from the B sublat- 
tice (o) to the A sublattice (•). This dimer configuration has 
winding numbers (0, 1) or (0, —1). 
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